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Abstract. Dynamics of quenching temperature is studied in pure and random Ising 
chains. Using the Kibble-Zurek argument, we obtain for the pure Ising model that the 
density of kinks after quenching decays as 1/ \fr with the quench rate of temperature 
1 /t for large r. For the random Ising model, we show that decay rates of the density of 
kinks and the residual energy are 1/ lnr and l/(lnr) 2 for large r respectively. Analytic 
results for the random Ising model are confirmed by the Monte-Carlo simulation. Our 
results reveal a clear difference between classical and quantum quenches in the random 
Ising chain. 



1. Introduction 

The change of a parameter across a phase boundary in a macroscopic system induces 
the dynamical phase transition of the system. If the changing speed of the parameter is 
sufficiently slow, the transition propagates over the whole system. However, as far as the 
changing speed is finite, symmetry breaking does not take place globally but does locally. 
It follows that spatial inhomogeneities emerge after the time evolution. Recent progress 
in experimental techniques enables to demonstrate such a dynamics across the phase 
transition and to observe imperfections in the ordered state. Greiner et al. pQ studied the 
dynamics across the quantum phase transition between the superfluidity and the Mott 
insulator in the optical lattice system. Sadler et al. [2] observed formation of defects after 
the quantum phase transition from the paramagnetic state to the ferromagnetic state 
in the atomic Bose-Einstein condensate. Weiler et al. [3] observed an evidence for the 
creation of vortices after the thermal phase transition of the Bose-Einstein condensation 
in the atomic Bose gas. 

The imperfection of the state after the evolution across the phase transition decays 
monotonically with decreasing the changing speed of the parameter. The decay rate 
depends on the choice of the parameter and the character of the associated phase 
transition. If the temperature is quenched in the classical system, the system undergoes 
the classical and thermal phase transition. If the strength of the quantum fluctuation 
is reduced at zero temperature, the quantum phase transition rules the character of the 
state. In the present paper, we study the dynamics near the thermal phase transition. 
Results are contrasted with those obtained for the quantum dynamics to reveal that 
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there exists a clear difference between the thermal phase transition and the quantum 
phase transition in the decay rate of an imperfection after the time evolution. Such 
quench dynamics across the phase transition are closely related to the dynamics of 
simulated annealing [I] and quantum annealing [5j [6J, [T] . It is an important issue whether 
quantum annealing performs better than simulated annealing or not [8]. 

The dynamics across the phase transition is well understood by the Kibble-Zurek 
mechanism [9], [10]. The scenario of the Kibble-Zurek mechanism is described as follows. 
Suppose that the system is driven from a disordered state to an ordered state by changing 
a parameter with a finite speed. The order parameter of the system in the initial state 
is uniformly zero. When the parameter comes close to the critical point, the healing 
time is so long that the parameter is changed further before the system attains the 
static state. Hence the system cannot evolves into the perfectly ordered state after 
the parameter passes through the critical point. The system after evolution consists of 
domains with different phases of the order parameter. The average of the domain size, 
i.e., the correlation length grows with decreasing the changing speed of the parameter. 
The growing rate is a universal function of the changing speed. 

The pure and random Ising chains are the models that permit us to investigate their 
dynamical properties analytically. Both of these models exhibit the quantum phase 
transition in the presence of the transverse field. The dynamics across the quantum 
phase transition in the pure Ising chain in the transverse field was first studied by Zurek 
et al. on the basis of the Kibble-Zurek argument [11] . They showed that the density p 
of kinks between ferromagnetic domains behaves as 

p Q (r)~l/v^ (1) 

with the quench rate 1/r of the transverse field. This result is confirmed by the analytic 
solution of the Schrodinger equation by Dziarmaga [12] • As for random systems, 
Dziarmaga applied the Kibble-Zurek argument to the quantum phase transition of 
the random Ising chain in the transverse field and obtained density of kinks decaying 
approximately as [13] 

pV)~l/(lnr) 2 , (2) 

for large r. Caneva et al. also derived the same decay rate [2], using the Landau-Zener 
formula and the distribution of excitation gaps at the critical point. Besides the density 
of kinks, they also estimated the decay rate of the residual energy. The result is given 
by 

EQJj) ~ l/(hir) C , C^3.4. (3) 

The study for the dynamics across the thermal phase transition in the Ising system 
is not necessarily sufficient. Laguna and Zurek have studied the Langevin dynamics of 
the order-parameter field in one spatial dimension [15]. However the model they studied 
does not correspond to the true Ising model in one dimension, because it involves an 
unphysical phase transition. Huse and Fisher have made a theory on residual energy 
after quenching temperature in classical random systems [16] . They regarded the system 
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with disorder as a collection of independent two-level systems, and derived residual 
energy which decays as 



Despite of Huse-Fisher's general theory on quenching in classical random systems, 
one cannot tell anything about comparison between dynamics across quantum phase 
transition and thermal phase transition. Comparison of eqs. ([3]) and (Hj) is obscure 
because of a lack of analytical support on eq. ([3]). Density of kinks tells nothing since 
it is not available in the classical dynamics. 

The results we obtain are summarized as follows. The density of kinks after 
quenching temperature in the pure Ising chain is given by 



for large r. As for residual energy, eq. (HJ) is reproduced for the random Ising chain. 
We emphasize that its derivation uses a manner different from the Huse-Fisher's theory. 
Our results reveal a clear difference in the decay rate of the density of kinks between 
the quantum quench and the classical quench. 

This paper is organized as follows. At first, we briefly review the Kibble-Zurek 
argument in the next section. After that, we study quench dynamics of the pure Ising 
chain and derive the decay rate of the density of kinks in sec. [31 We then reveal 
logarithmic decay rates of the density of kinks and the residual energy for the random 
Ising chain in sec. HI We also show results of Monte-Carlo simulation for the random 
Ising chain there. The paper is concluded in sec. [5j 

2. Kibble-Zurek argument 

Let us consider a ferromagnetic system with the critical temperature T c . In the Kibble- 
Zurek argument, the correlation length and the relaxation time of the system with a fixed 
temperature are quantities of importance. Both quantities are functions of temperature 
and increase with decreasing temperature toward T c . We denote the correlation length 
and the relaxation time by £(T) and T r (T) respectively. 

Now we consider quenching temperature with time t as T(t) — T c (1 — t/r) where 
time is assumed to evolve from — oo to r and 1/r stands for the quench rate. We 
assume that the system is in its equilibrium state initially. When the temperature is 
sufficiently high, the system almost maintains its equilibrium since the relaxation time 
is short. However, when the temperature is close to T c , the temperature decreases 
further before the system attains the equilibrium. Thus the system cannot possess the 
complete ferromagnetic order and contains domain walls when the temperature goes 
below T c . Once the domain structure forms, it should preserve until the temperature 
reaches absolute zero. The size of the domain is represented by the correlation length 



£ r C csO-)~l/(lnr) 2 . 



(4) 




(5) 



for large r. In the random Ising chain, it behaves as 
p c (r) ~ 1/lnr 



(6) 
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£ of the state when the temperature passes T c . An argument described below provides 
an estimation of £ for a given r. 
We introduce an equality: 

r r (T(t)) = \t\. (7) 

This equality defines the time t at which the relaxation time is equal to the remaining 
time to the critical temperature. At a later time until t — r, the system cannot attain 
the equilibrium since the relaxation time is longer than the remaining time. Suppose 
here that the system stays in the equilibrium at T = T(t) and does not evolve any more 
after t passes t. Then the correlation length of the state at t = r is approximated by 
£ ~ £(T). Since one can express T in terms of £ from the expression of £(T), the left 
hand side of eq. ([Zj) is written in terms of £. The right hand side, on the other hand, 
is written as r\T — T c \/T c , which is also expressed in terms of £. Thus we obtain the 
equation of £ from eq. ((Tj). The solution of this equation yields £ as a function of r. 



3. Pure Ising chain 

We consider the simple pure Ising model in one dimension: H = ffi&i+i- Although 
this model does not exhibit the phase transition at any finite temperature, the ground 
state possesses the complete ferromagnetic order. Hence one can regard the critical 
temperature as T c = 0. Denoting the inverse of temperature by (3, an expression of the 
correlation length is given by 

£(T) = l/lncoth/3 « ^e 2/3 , (8) 

where the lattice constant is assumed to be the unit of length. We note that the right 
hand side is the expression valid at low temperature, i.e., Tel. In order to discuss the 
dynamics of the present system, we assume the Glauber model [IT]- Then, the relaxation 
time for a fixed temperature is given by 

r r (T) = 1/(1 - tanh2/3) « « \t(T)\ (9) 

where the approximation signs are valid at low temperature. Thus the correlation length 
and the relaxation time grow with decreasing temperature toward T c = 0. 

From now we discuss the dynamics of quenching temperature according to the 
Kibble-Zurek argument. We assume the quench schedule: 

T(t) = -t/r (10) 

instead of the one in the previous section because of T c = 0. We also assume that 
the time t evolves from — oo to and the inversed quench rate r is large, i.e., r ^> 1. 
Equation ([7j) defines the approximate time t at which the evolution of the system stops. 
Using eqs. (jSJ) and (|10p . the time t is written as t « — It j In 2£, where £ is the correlation 
length at T = T(t). We remark that the low temperature expression of £ is allowed as 
far as r is large enough because T is small. Equations ((Tj) and (jHJ) yield an equation 
of £ as £ 2 ~ 4t/| ln2£|. This equation cannot be solved analytically. However ln£ is a 
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gentle function of £ compared to £ 2 . Hence £ is almost proportional to y/r. The inverse 
of correlation length corresponds to the density of kinks approximately. It follows that 
the density of kinks in the final state is estimated as 

1 (I ln2el) 1/2 , N 

p^-^ = ^ ^ — . 11 

Thus one finds that the density of kinks of the final state is proportional to 1/ \fr as far as 
the logarithmic correction is ignored. The logarithmic correction is not essential indeed. 
To verify this, we consider a modification of the Kibble-Zurek argument as follows. We 
may employ another equality, r r (T(t)) = r, instead of eq. (J7J). This equality defines the 
time t at which the relaxation time exceeds the inverse of quench rate, we can consider 
that the evolution of the system almost stops at t and construct the argument same 
as that in the previous section. By this argument, we obtain p ~ without the 

logarithmic correction. 

4. Random Ising chain 

The random Ising chain is represented by 

H = -J2Jm<7i+i- (12) 

i 

In our study, the coupling constant {J{\ is drawn randomly from the uniform distribution 
between and 1, namely P(Ji) = 1 for Jj G [0, 1] and P{Ji) = otherwise. This model 
corresponds to the one studied in refs. [13l [H] . 

The correlation function between sites i and % + k with fixed { Jj} in the equilibrium 
at a fixed temperature is given by (crjCTj+fc) = n}t=i _1 tanh. (3 Jj. Taking the average 
over randomness, the correlation function in the thermodynamic limit is obtained as 
[(o~iai + k)]a,v — (lncosh/3//3) fc . From this formula of the correlation function, one can 
obtain an explicit expression of the correlation length: 

£(T) = [ln(/5/ In cosh w (3/ In 2. (13) 

Note that the right hand side is the low temperature expression. 

The energy of the system with fixed { Jj} is written as (7i) = — J2i -h tanh /3Jj. The 
average over randomness yields an expression of the energy per spin at low temperature 
in the thermodynamic limit: 

e = hm — — - — « h tht — • (14) 

Jv^oo N 2 (3 2 24 v ; 

We remark that the ground state energy is — |. 

The relaxation time is available in ref.[18] by Dhar and Barma. It is given by 
Tr = 1/(1 — tanh 2/5). The low temperature expression is 

r r (T) ~ l -e^ ~ I e ( 41n2 ^ T ). (15) 

As is the case with the pure Ising chain, the critical temperature of the present model 
is T c = 0. 
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Let us consider that the temperature is lowered according to the schedule given 
by eq. (ITUl) . We impose eq. ([7j) to define the time i at which the evolution keeping 
equilibrium breaks. From eq. ( |T3l) . the time relates with the correlation length by 
\t\ = r/(£ln2). Applying this relation and eq. (fT5|) to eq. ([7]), we obtain an equation 
off: 

i= J-flnr-lni^V (16) 
^ 4 In 2 V 2 J K ' 

This equation cannot be solved analytically. However, since (ln£)/£ — > for £ — > 0, 
we find that £ is almost proportional to lnr when r 3> 1. Equation f[TB"j) leads to an 
estimation of the density of kinks, 

1 41n2 , , 

n~^~ s . 17 

£ lnr-ln^ 

The second term in the denominator is negligible for a sufficiently long r as mentioned 
above. Hence eq. (JBj) is derived. 

The residual energy is also estimated from the energy at T = T. Using eqs. ( fl5l) 
and (ITUl) . eq. (I7j) is rewritten as 

\^ = t/$, (18) 

where we defined (3 = 1/T. This equation is followed by (3 = |lnr — |ln(/3/2). 
Substituting this for (3 in eq. (TT41) . we obtain the residual energy per spin as 

2tt 2 1 



.. • (19) 

3 (lnr-ln(/3/2)) 2 V V 

Since the second term in the denominator is negligible for large r, hence we obtain eq. 

(©■ 

We next consider a logarithmic schedule: 

T(£) = T ° T - , (20) 

W l + aln(-^f)' V 7 

where T and a are positive numbers. In this schedule, the temperature is reduced from 
T at t = —T t to at t — 0. Using eq. ( |20|) with eqs. (JTj) and (TT5I) . one obtains 



the equation of f3 as ^e 4/3 = T rexp{l/a — (T /a)/3}. This equation can be solved 
analytically and yields $ = ln(2e 1/a T r)/(4 + T /a). From eq. dTJ), one obtains the 
expression for density of kinks as 

1 (4+ 2>) i n 2 

P ^lnr + ln(2T ) + i- (21) 

This expression is reduced to eq. (jSJ) for r — > oo. The expression of the residual energy 
per spin is obtained as 
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Figure 1. The density of kinks after cooling with the linear schedule. Square symbols 
are obtained by the Monte-Carlo simulation. The fitting curve is given by p = q/£ 
with the correlation length £ determined by eq. (|16p with r = pr mc . Parameters p 
and q are determined by denotes of the method of least squares. 

which yields eq. fl4]) for r — > oo. These results imply that asymptotic behaviors of the 
density of kinks and the residual energy for r — ► oo are insensitive to the schedule of 
quenching temperature. 

If the distribution of J^ has a finite positive lower-bound, decay rates of the density 
of kinks and the residual energy are the same and obey the power law. To show this, 
we suppose G [Jo, 1] with < J < 1. Then the correlation length is written 
as £ ~ f3AJe 2 P Jo , where we defined A J = 1 — Jo. The energy per spin is given by 
e w e g + pXjJo e ~ 2 ^ J ° 5 where e 9 = — (1 — Jq) is the ground-state energy. The 
relaxation time does not change by introducing Jo [18]. It follows that the condition 
in the Kibble-Zurek argument, eq. (j7j), yields the same equation of (3 as eq. (ITS]) . 
Using its solution for large r, We obtain e ms ~ Jop and p ~ r _Jo//2 /lnr. As for the 
quantum quench, the introduction of a finite positive Jo does not change the universality 
of quantum phase transition [19]. It follows that the decay rate of the density of kinks 
is logarithmic and given by eq. ([2]). Hence the classical quench reduces the density of 
kinks faster than the quantum quench in this case. 

We confirm results of the random Ising chain by the Monte-Carlo simulation for 
systems with 500 spins. The temperature is lowered according to the linear schedule. 
We choose the initial condition for the temperature as T = 5 at t = — 5r. The coupling 
constant J^ is drawn from [0, 1] uniformly. In order to take an average with respect to 
randomness of the system, we generated 100 configurations of coupling constants { J{\. 
For each configuration, we perform quenching temperature 500 times. 

Square symbols in Fig. [JJ and Fig. [2] show the density of kinks and the residual 
energy per spin respectively obtained by the Monte-Carlo simulation. The density of 
kinks is defined by p = [j^ I]j(l — (c r ?c r j+i))]av ) where (■ • •) denotes the expectation value 
with respect to the state after simulated annealing and [• • -] av means the average over 
configurations of coupling constants. 

In order to obtain fitting curves for the density of kinks and the residual energy per 
spin, we need to modify eqs. (fTTI) and (fT9l) . First, we have to care about the difference 
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Figure 2. The residual energy obtained by Monte-Carlo simulation (square symbols). 
The fitting function is given by e rcs — ^/(rP) 2 with 0, the solution of eq. (fT5)) with 
t = pT mc . The fitting parameter r is determined by the method of least squares with 
p given from fitting of density of kinks. 

in the unit of time between the Glauber's dynamics and the Monte-Carlo dynamics. 
Then we bring up the relation, r = pr mc , between the inverse of quench rate r in the 
Glauber's dynamics and r mc in the Monte-Carlo dynamics, where p is an adjustable 
parameter. Next, we relate the density of kinks p to the correlation length £ by p = g/£, 
where £ is the solution of eq. (fTBl) . The parameter q tunes the inverse of correlation 
length to the density of kinks. Finally, we propose an ansatz that the residual energy is 
represented by e Tes = (tc 2 /24)/(r/3) 2 with an adjustable parameter r, where $ is given 
from eq. (Fl8|) with r = pr mc . 

Parameters p and q are determined by the method of least squares described as 
follows. The Monte-Carlo simulation yields a set of data (pi, t™ c ), where pi is the mean 
value produced by Monte-Carlo simulations. The fitting function given by eq. (fP7j) but 
p = q/£ and r = pr mc does not yield p for given r mc analytically. Then we regard 
r mc as a function of p and define the error between Monte-Carlo data and the fitting 
function by S(p, q) = Z)i(l nrmc (Pi) ~^ aT i Clc ) 2 / (T i where r mc (p;) is the value of the fitting 
function and r™ is given from the Monte-Carlo data for pi. of is the dispersion of p«. 
We assume that the relative values of the dispersion between different i's are the same 
in p and In r mc . By minimizing S(p, q) , we fix the value of p and q. The errors of p and q 
are given by a 2 = -^ttitt-ft and a 2 = -^rCf-FT, where m is the number of Monte-Carlo 

° J P m—2 aq z D 1 in— 2 ap z D ' 

data and D = ^4|r4 — (-tNt) 2 . The other parameter r is determined by the method 
of least squares with the error S'(r) = Z)i( £: (' r i mc ) — e i) 2 ' between the value of fitting 
function, eq. (TT9j) . with /3j given by eq. ([18]) with = pr™ c and the value of Monte- 
Carlo data £j. o 2 ei is the dispersion of £j. p is fixed by fitting of (p, r mc ). The error of r 
is given by of = ^^-(^i-) -1 - The obtained values of p, q and r are g ~ 0.241 ± 0.001, 
p w 22.0 ± 0.9, and r w 2.14 ± 0.16. 

Figure [1] and [2] shows that results of Monte-Carlo simulation on the density of kinks 
and the residual energy are excellently fitted by the curves made from eq. (TP7I) and (TT9l 
respectively. Therefore the analytic results on the basis of the Kibble-Zurek argument 
is confirmed by the Monte-Carlo simulation. 
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5. Conclusion 

We studied the dynamics of quenching temperature of pure and random Ising chains, 
on the basis of the Kibble-Zurek argument. We showed for the pure Ising chain that the 
density of kinks after quenching decays as 1/ -Jr for large r. As for the random Ising 
chain with Jj e [0, 1], the density of kinks and the residual energy decay as 1/ lnr and 
l/(lnr) 2 for large r respectively. Results for the random Ising chain were confirmed by 
the Monte-Carlo simulation. Comparing our results on the density of kinks with known 
results for the quantum quench, densities of kinks after the classical quench and the 
quantum quench decay with the same power of r in the pure system. As for the random 
Ising chain, the power of In r by the quantum quench is twice as large as that by the 
classical quench. The difference between the quantum quench and the classical quench 
is substantial. 

The classical quench and the quantum quench toward the ground state correspond 
to simulated annealing and quantum annealing respectively of an optimization problem. 
The random Ising chain studied in the present paper provides an optimization problem 
with the trivial solution. However it is a non-trivial problem for simulated annealing 
and quantum annealing since their dynamics toward the solution respond to randomness 
and exhibit the slow relaxation. Our results reveal that the random Ising chain is a 
solid example for which quantum annealing certainly performs better than simulated 
annealing. 

The author thanks T. Caneva, G. E. Santoro, and H. Nishimori for fruitful 
discussions. The present work was partially supported by CREST, JST. 
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